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The critical temperature T c of an interacting Bose gas trapped in a general power-law potential 
V(x) — Ui\xi\ Pi is calculated with the help of variational perturbation theory. It is shown that 
the interaction-induced shift in T c fulfills the relation (T c - T°)/T° = Di(n)a + D'(n)a 271 + 0{a 2 ) 
with T<? the critical temperature of the trapped ideal gas, a the s-wave scattering length divided 
by the thermal wavelength at T c , and r\ — 1/2 + the potential-shape parameter. The terms 

D\(rj)a and D'(rj)a 2ri describe the leading-order perturbative and nonperturbative contributions to 
the critical temperature, respectively. This result quantitatively shows how an increasingly inho- 
mogeneous potential suppresses the influence of critical fluctuations. The appearance of the a 2ri 
contribution is qualitatively explained in terms of the Ginzburg criterion. 

Q . PACS numbers: 03.75.Hh,05.30. Jp,64.60.Ak 

I. INTRODUCTION 

The realization of Bose-Einstein condensation (BEC) in dilute atomic vapors has renewed interest in the critical 
properties of weakly interacting Bose gases and, in particular, their transition temperature T c . Important recent 
work in this area concerns the role of the external trapping potential. For the homogeneous Bose gas, the shift in 
T c caused by s-wave contact interactions is in leading order completely due to long-wavelength, critical fluctuations 
that have to be described nonperturbatively. It is now established that these fluctuations lead to a linear increase of 
the critical temperature with the s-wave scattering length a, if the particle density is fixed 0, U 03, IH . A harmonic 
trapping potential, on the other hand, suppresses the critical long-wavelength fluctuations and reduces the fraction of 
£j ' atoms taking part in nonperturbative physics at the transition point. As a result, the leading-order shift in T c can be 
calculated by simple pertubative methods, for instance the mean-field (MF) approximation to the Landau-Pitaevskii 
equation 0,0]. Interestingly, the shift here is a decreasing linear function of the scattering length a. 

Work on the critical temperature of Bose gases has so far been mainly concerned with homogeneous and harmonically 
i trapped systems. As outlined above, in these two situations very different physical mechanisms determine the shift of 
T c . This observation naturally motivates an investigation of the crossover between the two cases. In this paper we shall 
t-H , interpolate between these limits by studying the condensation in a general power-law potential, whose parameters can 
be varied continuously. In this way we shall obtain a deeper understanding of how the increasing inhomogeneity of the 
potential suppresses the critical fluctuations and changes nonperturbative into perturbative physics. The power-law 
potentials under study are given by 



(1) 



with Ei and Li denoting energy and length scales. If all powers pi are set equal to 2, we recover the harmonic potential, 
whereas in the limit of all pi diverging, T^(x) approaches a box shape characteristic of the homogeneous Bose gas. 

First investigations of the crossover behavior of the critical temperature in these potentials were recently reported in 
Refs. H@. In 0, the shift in T c for these potentials was determined within mean-field theory in the thermodynamic 
limit. Extending earlier first-order calculations [9lll0|. it was shown that up to second order in the scattering length 
\ a, the MF shift has an expansion 

^-^ = D[ MF \ V )a + D[ MF) ( V )a^ + + o(a 2 ) (2) 

with T° the critical temperature of the noninteracting gas, and a = a/ '\t c the scattering length measured in units of 



the thermal wavelength At = y 2irh 2 /ksTm for particles of mass m at a temperature T. The exponent of the second 
term is twice the shape parameter of the potential, 
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so that 2ij grows from 1 to 4 as the shape changes from homogeneous to harmonic. The coefficients D'i MF \ [t]) and 

D^ F ' (rf) are respectively given explicitly or through simple quadratures. These results provided first insights into 
the crossover in the behavior of T c between homogeneous and inhomogeneous potentials. However, since mean-field 
theory does not account for critical fluctuations, it can only provide a rough first estimate, especially in the quasi- 
homogeneous regime, and needs to be improved by more sophisticated methods. 

One possible pathway for taking critical fluctuations into account was explored in Ref . @ , where the shift in T c was 
calculated with the help of a renormalization group (RG) method initially developed for studying the harmonically 
trapped gas pdT |. The results were found to be in good qualitative agreement with mean- field theory. As a main 
advantage, the RG approach employed in that work gave a simple and transparent tool to compute the critical 
temperature for a wide range of potential shapes and interaction strengths. A disadvantage was, however, that the 
results were mainly numerical and rendered only limited physical insight into the system. Furthermore, the calculation 
required several unsystematic approximations which are difficult to improve. 

The purpose of this paper is to present a more systematic approach to the problem by making use of field-theoretic 
variational perturbation theory (VPT). VPT is a powerful resummation method for divergent perturbation series 
[To| , which has been extended to quantum field theory and its anomalous dimensions in Ref. 0, ^| . It has led to 
a prediction of the critical exponents of superfluid helium with un prec edented accuracy |l3l Il4| , as confirmed for the 
exponent a of the specific heat of helium by satellite experiments |15| . 

In the context of BEC, field-theoretic VPT has been applied successfully to determine the shift of the critical 
temperature of the homogeneous Bose gas from a five-loop perturbation expansion |T3 |. extended recently to six 
and seven loops in |17| |. In the present work we shall describe the trapped Bose gas in the thermodynamic limit, 
in which the trap is so wide that we may apply, as in the local-density approximation (LDA). The system is 

treated as locally homogeneous at any point x with an effective chemical potential /i ff = fi — V(x), where /i denotes 
the global chemical potential. In this way, we can make contact with the high-order perturbative loop expansions 
that were derived in Refs. 0, 0] for classical three-dimensional (/> 4 -theories of homogeneous systems. Because of 
dimensional reduction 0] , the effective classical theory 0, 0] can be directly used to describe critical properties of 
the (quantum-mechanical) Bose gas below second order in the scattering length 0, • 

In this work we shall combine the high-order loop expansions with the LDA to derive a perturbative expansion for 
the particle number of the trapped system in powers of a = a/Xr- From this expansion we extract, with the help 
of field-theoretic VPT, the critical particle number and the shift of T c . The main results are: (i) For small a, the 
shift of T c is shown to exactly follow a behavior (T c — T°)/T® = Di(rj)a + D'(r])a 2v + 0(d 2 ) in generalization of the 
mean- field result @. The term proportional to a 2v represents the leading nonperturbative effects, whereas Di(rf) can 
be calculated perturbatively as discussed in 0,0. The second-order contribution, which we will not study in detail, 
contains terms proportional to a 2 and a 2 In a. (ii) We compute the coefficient D'(rj) for r\ < 1 using VPT, and in 
this way arrive at a quantitative description of the behavior of T c below second order in the scattering length, (iii) 
Following , we give a qualitative explanation for the a 2ri behavior of T c based on the Ginzburg criterion [5(J . 

The paper is organized as follows. In Sec. II, the physics of ideal Bose gases in power-law potentials is briefly 
reviewed. In Sec. Ill we show with the help of general scaling arguments why the critical temperature obeys a law of 
the form and give a physical interpretation for the appearance of the nonanalytic term. Section IV contains the 
calculation of the nonperturbative coefficient D'(rj), and illustrates the behavior of the critical temperature by means 
of a numerical example. Summary and conclusions are given in Sec. V. 



In this section we summarize some known properties of ideal Bose gases in power-law traps that are needed later. 
The notation follows Refs. Q and [2l|. We consider a system of N ideal bosons of mass m trapped in the power-law 
potential characterized by the shape parameter 77 introduced in Eq. (J3J). We define the characteristic volume by 



where I(pi) = T(l/pi)/pi and T(z) denotes the usual gamma function. The quantity Vchar provides an estimate for 
the volume occupied by the one-particle ground state in the trap. 

For calculations in the local-density approximation it is useful to convert spatial integrations involving the trap 
potential into energy integrations according to the rule J d 3 x f[V(x)] — J dep(e)f(e). In this way, we can easily deal 
with power-law potentials of any type. The density /5(e) is the area of the equipotential surface V'(x) = e. As shown 



II. IDEAL BOSE GASES IN POWER-LAW TRAPS 
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p(s) is given by 

v 2( V +l)/3 / x 7,-1/2 

As announced above, we shall treat the thermodynamic limit, in which N and V^har go to infinity at fixed 
^/Khar +1 ^ 3 - The equation of state for the ideal Bose gas above the condensation point is then given by p3. l23| 

w =(^(S) ,+, ^ +w V.w. w 

where /3 = 1/ksT denotes the inverse temperature, \x the chemical potential, and z — exp(0fi) the fugacity. The 
Bose-Einstein functions Ca(^) = EfcLi z k /k x are polylogarithmic functions [2J|. The spatial density distribution of 
the gas is determined by the relation 

n(x) = A T 3 C 3/2 (e^-^)l) . (7) 

The condition for Bose-Einstein condensation N — J d 3 x n(x, n = 0) can be obtained from Eq. © by setting z = 1 
or directly from Eq. J7J and reads [2lJ 

where we have replaced by Riemann's ^-function £(,,) = Y^Li ') an( i = l/feTc denotes the inverse 

critical temperature of the ideal gas. 



III. GENERAL BEHAVIOR OF CRITICAL TEMPERATURE 



We now turn to the interacting Bose gas where we assume, as usual, that the interaction is effectively a delta-function 
potential which is completely characterized by the s-wave scattering length a. In the local-density approximation to 
the thermodynamic limit, the trapped particle number N at given temperature T and chemical potential \x can be 
calculated from the integral 

N{T,n) = J d 3 «tr(x;/i,T)^ J d 3 x n hom (^i — V(x), T), (9) 

where n tr is the trapped particle density and nh om the density of the homogeneous gas. As we briefly explain at the 
end of this section, we expect the LDA to be applicable if the condition A 2 ^ /a <C is fulfilled, where At c denotes 
the thermal wavelength at the condensation point and V^ar is defined in Eq. Q). Obviously, for a fixed At c this 
condition can always be met by making the trap wide enough (and increasing the particle number accordingly) . 

In the following, we want to apply Eq. to explain why the critical temperature follows a behavior as indicated 
in Sec. I. A more detailed calculation will be presented in Sec. IV. Consider the perturbation expansion of the 
homogeneous density nh om in powers of a = a/ At 

n hom (T, A/x) = n^jT, A/i) + n^jT, A^)a + n^jT, A^)a 2 + O (a 3 ) , (10) 

where A/i = /i c — \i is the negative distance of the chemical potential /x from its critical value /x c at temperature T. 
Our definition of Afj, ensures that it is positive above the transition. 

The omitted higher-order terms in the expansion (|10fl depend on the details of the particle interactions. Being 
interested in contributions below second order, we can disregard these details and work only with a contact interaction. 
Note that in Eq. (HU| we have neglected logarithmic terms appearing at second and higher order in a. These terms 
enter via the critical chemical potential \i c which contains a contribution proportional to a 2 In a |(| . As the calculation 
of Sec. IV shows, the omitted logarithmic terms are not relevant for determining the shift in T c below second order. 

From Eq. (|1U|) . we can convince ourselves that there must exist a perturbative second-order contribution to the 
critical particle number. This justifies the inclusion of a term proportional to a 2 in the general expression for the 
shift in T c , as mentioned at the end of Sec. I. Indeed, since the expansion (|10|l becomes arbitrarily accurate when we 
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go sufficiently far away from the critical region, we can split the spatial integral in Eq. @ into a part near the trap 
center and a remainder: 

N{T,fi c )= [ d 3 xn hom (^ c -V(^),T)+ [ d 3 xn hom {fi c -V(x),T). (11) 

JV(x)<Vo JV(x)>V 

Here, Vq denotes an energy above which the perturbative expansion of the density becomes accurate. Inserting the 
expansion (|ll)f) into the second integral, we see that N(T,/j, c ) contains a contribution which is of second order in the 
scattering length. Since the expansion for nh om contains a term proportional to a 2 In a as mentioned above, there will 
also be such a contribution in the exact second-order result (compare to Ref. Q for the harmonic case). 

The expansion coefficients and nP^ in the perturbation series (|10|l remain finite in the critical limit Afi — > 0. 
From Eq. 10 we thus obtain well-defined perturbative contributions to the critical particle number in zeroth and first 
order in a. The zeroth-order contribution is just the critical particle number of the ideal gas. However, all other 
coefficients ^iom' * — 2, are infrared-divergent in the critical limit of A/i — > where they behave like 1/y/Afi . 
Nevertheless, as shown in Refs. [Til Il7j. we can make use of resummation techniques to extract information about 
critical properties from the expansion. If we focus on effects below second order in the scattering length, it is sufficient 
to consider only the leading divergence ftj^ at each order, i.e., = + o((A/i)( <_1 )/ 2 ). This leads to the 

power series 

oo oo , „ s i— 1 

A<UT, A M ) ee £ (r, Atf* = a £ bi (-^=) . (12) 

Above the transition where the chemical potential [i is smaller than /i c so that A/x > 0, we insert <|12[) into © and 
obtain the leading-order divergent contribution to the trapped particle number 

AN div (T, Am) = / rf 3 a:A<;(r, A M + V(x)) =a j d 3 x £ b t ( = ) . (13) 



Converting the spatial integral into an energy integral with the help of JHJ and performing this integral with analytic 
regularization, we obtain 



/oo 
dee r '- 3/2 a^2b i 

OO 



i=2 



y/3(AM + e)/ 
a Y T(i/2-r)) 



jpEil) r(»72-i/2) 



= ^ Vh \vm)' (14) 

where irrelevant constants have been absorbed into the coefficients A and C. Note that the factors T(i/2 — 77) in the 
coefficients cause divergences for r\ — > 1. We ignore this issue for the moment and defer its discussion to Sec. IV. 
The main property of (|14|) is that the function h\ in the final expression depends only on the ratio a/y//3Afi. If the 
number of particles is to remain finite in the critical limit A/x — > 0, the limiting behavior of hi(d/\//3Afi) must be 
hi(a/y/PA~}l) cx (a/V/3A7I) 2 '>. It follows that the critical AN div behaves like a 2 ''. 

Combining this result with the perturbative first-order contribution mentioned above, we find that the change 
AiV cr it in the critical particle number at fixed temperature is given by 

A ^ = C 1 (7 ? )a + C'(r ? )a 2 "+0(a 2 ) (15) 

- /V crit 

with AT? rit the critical particle number of the ideal gas, and C\{rj), C'(rf) proportionality constants depending on the 
potential shape. From this behavior we immediately deduce the change of the critical temperature at fixed particle 
number to behave like 

AT 

-^= Dl (r ] )a + D'(r 1 )a^ + 0(a 2 ) (16) 



with coefficients D\{rj) and D'(rj) which follow trivially from Ci(ry) and C'irj) [compare with Eq. (|25|l below]. 
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FIG. 1: Schematic diagram of shift in T c for a fixed small scattering length a as a function of the potential shape parameter 
r\. The upper curve shows the full result below second order according to Eq. 1161 . the dashed lower curve only displays the 
perturbative contribution linear in a. The nonperturbative contribution decreases fast with growing inhomogeneity. 



To discuss the physical contents of Eq. we anticipate some of the results of the next section, and schematically 
show in Fig. ^ the behavior of the critical temperature at a fixed, small value of a as a function of the potential shape 
parameter r\. The full curve shows the combined perturbative and nonperturbative contributions, whereas the dashed 
curve displays only the perturbative (linear) term. In the homogeneous limit we have 2rj = 1, so that the shift below 
second order in a is purely linear as we expect from earlier studies 0,0, H^|- In this case, both contributions [i.e., 
Di{rj)a and D'(r])a 2ri ] are of comparable size at any value of a. 

The situation changes when we enter the inhomogeneous regime where 2r) > 1. As displayed in Fig.^ for sufficiently 
small, fixed a and growing rj the nonperturbative contribution D'(rf)a 11 rapidly becomes very small compared to the 
perturbative term. This kind of behavior is independent of the detailed form of D'(j)) [note that in Fig.^ we ignore 
the (unphysical) divergence of our approximation l|ltj|) in a narrow vicinity of 77 = 1; as discussed in Sec. IV, this 
is expected to be remedied in a higher-order expansion]. Equation (|lbl describes quantitatively how the growing 
inhomogeneity of the potential reduces the influence of critical fluctuations on the transition temperature. 

Following the arguments of Ref . [|| , we can also give a physical explanation for the appearance of the a 2ri term by 
estimating the fraction of atoms actually taking part in nonperturbative effects. For simplicity, consider the potential 
V(x) = Vo|xi/ro| 7 . Treating the Bose gas above the transition point as classical, we find the mean-square width 

^2x_-2 r(3/7) , 1 = 1,2,3, (17) 



x 11 u (/3Vb) 2 /7r(l/ 7 ) 
i.e., the cloud radius behaves like 

r-doud- (/JVb)1/7 . U«j 

From the Ginzburg criterion [2£j it follows that nonperturbative effects only arise at (local) chemical potentials fi e g, 
for which 

h 2 a 2 , , 

Mc-^off < — jr. (19) 

III Arp 

Invoking the local density approximation, this means that the nonperturbative region around the trap center has a 
radius of about 

n a x 



■ (20) 



The fraction of atoms within this nonperturbative spatial region is given by 



np 



3 / „ \6/ 7 



r c ioud / \ At / 



(21) 



However, not all atoms within this region actually take part in nonperturbative physics 0- From the homogeneous 
system we infer that only a fraction a/ At are actually involved in these effects. For the trap, this means that the 
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fraction of "nonperturbative" atoms scales like 



a 



np 



(6/7) + l / „ \ 2-) 

(22) 



aV 6/7)+1 - (— ^ 



At V r cloud/ \At/ \At/ 

which explains the appearance of the nonanalytic term in Eq. i|15[l • 

At this point, it is also convenient to explain the estimate for the validity of the LDA given above, again following 
the arguments of Ref. ||. Nonperturbative effects involve fluctuations with wavelengths A^/a and larger. For the 
LDA to be applicable, the size r np of the nonperturbative region around the trap center should thus be much larger 
than Afi/a. Since Eq. (|2Ul can be rewritten as 



At 



2/7 

F (2/ 7 +l)/3 
char 



by using Eq. this condition immediately implies that Ay/a -C K^.. In other words, the extension of the ground 
state has to be much larger than the minimum length scale A^/a for critical fluctuations. With the help of Eq. JSJ, 
we can also rephrase this statement as follows. With condensation taking place at temperature T, the trap has to be 
sufficiently wide, so that the critical atom number fulfills the condition 

aN l/{2r l +2) A 



IV. CALCULATION OF CRITICAL TEMPERATURE 



After the general discussion of Sec. Ill, we now turn to the actual calculation of the nonperturbative coefficient 
D'{rj). Our result will be approximate in two ways: first, we use some approximations to calculate the coefficients of 
the weak-coupling expansion for the trapped particle number. Second, the number of terms is limited to seven, and 
the evaluation via VPT leaves an error. However, due to the stability and the typically exponentially fast convergence 
of VPT, our results should provide a satisfactory representation of the true behavior. 

To simplify the notation for the following calculations, we introduce the reduced homogenous density function 

/(/3(/x c - n)) = A T n hom (/3/i) = A T n hom (/3/i c - (3(fj, c - //)). (23) 
This quantity implicitly depends, of course, on the reduced scattering length a = a/Xr, as indicated in the expansion 

Under the assumption of LDA, f(A + (3V(x)) with A = /3A/i > equals the reduced particle density of the trapped 
Bose gas at point x. For A = 0, it describes the critical density profile in the trap. We now define the integral 

/•oc 

I V (A) = dvv^^fiA + v). (24) 



Using Eqs. © and this can be recognized as a rescaled version of Eq. ©. The shift in the critical temperature 
for a fixed particle number is then given by 0, @ 

r -V +, _ Ch+D (25) 



t<o>J /,(o)/r( 1) -i/2)- 

As sketched in the previous section, we shall first derive a perturbative expansion for f([3(/i c — //)) in powers of a. For 
A — > 0, the zeroth- and first-order terms of this expansion remain finite, whereas the higher-orders terms suffer from 
infrared divergences. The zeroth and first order can thus be directly inserted into Eq. (|24[) and their contribution read 
off at A = 0. For the higher-order terms, we focus only on the leading-order divergence. This leads to an expansion in 
terms of a/y/A + v. Inserting this result into l|24f) and performing the integration over v, we obtain the weak-coupling 
expansion for I V (A) [compare with Eq. I|14(l ]. This expansion is finally resummed to find the coefficient D'(rf). 

Let us temporarily return to the more familiar unsealed quantities to outline further details of the calculation. 
Using the LDA and Eq. (J5J), the trapped particle number is given by 

•** ^hom 

(7>-V(x)) 

= (^)" 1/2 r ^^Kir l)/ 7 de ^ 3/2 ^( T ^- £ )- ( 26 ) 
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FIG. 2: Diagrams contributing to iV' 2 ' up to order a 4 . The crosses denote the joining of two Green functions. 



We now express nh om (T, fj, — e) in terms of the Green function of the interacting homogeneous system: 

1 

G h om(k,w„;/x,T) = — 



<>n - [£k- M+ 7i£(k,w„;At,T)]/7i' 

where £k = ?i 2 k 2 /2m with momentum k, uo n — 2im/f3Ti with integer n are the Matsubara frequencies, and 
E(k, uj n \ /i, T) is the proper self energy. This leads to 



N(T,v) = - -3 



RV r(?7 - 1/2)^(2tt) 3 

x / dee n ~ 3/2l Sf - (27) 

J ^Tkiun- (e k -Ai + e)/R-E(k,w n ;^-e,T)' 

The symbol ^J fc with fc = (k, w n ) denotes integration and summation over all momenta and Matsubara frequencies. 
The last term in (|27[) is conveniently rearranged to 



& 



[£k - (M - e) + ftE(0, w„; ^ - s, T)]/h - [E(k, w„; n - e, T) - E(0, w n ; M - e, T)] ' 
In this expression, we use the "mass-renormalized" Green function 

1 



(28) 



G (k,cj n ;fj,,T) 



iu n - [e k - fx + 7iE(0, w„; ^, T)]/7i 



as the free Green function for a perturbative expansion. In this way, we obtain two different contributions to the 
trapped particle number N(T, (i) = A^ (1) (T, (i) + 7V( 2 )(T, /x), namely, 

(i) the zero-order term 

NW(T,n) = C [ dee^/zyf — L- — - (29) 

V ; J 4^ n -[£ k -(/i-e) + ftE(0,w n ;/i-£,T)]/7i v ; 

with C the prefactor in front of the integral of Eq. (|27J) and 

(ii) higher-order contributions pictured by the Feynman diagrams in Fig. up to five loops. Using the zero Mat- 
subara frequency contributions of these diagrams, we will calculate a perturbation series for the leading infrared- 
divergent contribution n^^ v ^ (T, /i) to the homogeneous density. From this we obtain the second contribution 

]Sf( 2 ' dlv )(T,fx) = J d 3 x n[ 2 ^ lv '(T, ^ — U(x)) to the trapped particle number. We note that the nonzero Matsubara 
frequency modes are expected to contribute only in second order in a to the critical particle number 0. 

First we discuss AT' 1 ) (T, fi). The relevant diagrams for the perturbative evaluation of Eq. I|29ll are shown in Fig. 
03 up to three loops. As in Fig. [2 an n-loop diagram contributes to order a™ -1 to the perturbation series. The 
contributions of zeroth and first order in a are very easily calculated. Since they are convergent in the limit of A — ► 0, 
we find up to first order in a: 



1,(0) = rfa-l/2)Cfa + l)+o / dvv^ 3 /HF 1/2 (v)[C(3/2)-F 3/2 (v)}+... 

Jo 

= rfa - 1/2)C(7? + 1) + T(ry - l/2)Cifa)a + . . . (30) 
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|P=0 



FIG. 3: Diagrams contributing to N^'(T, fi) up to order a 2 . The sunset subdiagram is taken at the external momentum p = 0. 



This expression is equivalent to the first-order mean-field results of Refs. @- The perturbative contribution is 
nonzero in the homogeneous Unfit r? = 1/2 013. 

The higher-order diagrams are divergent in the limit of A — > 0. We shall not use the whole set of diagrams for 
our calculation, but restrict ourselves for simplicity to the Hartree-Fock (HF) approximation in which we take only 
diagrams into account that consist purely of simple bubbles (such as the first four in Fig. Unfortunately, it 
seems difficult to estimate the consequences of this approximation or to even go beyond it, but we expect that it 
captures the main features of the actual behavior. At any rate, the HF approximation is interesting in itself since 
the resummation can be done exactly and provides a nice illustration of our approach. We also find that it leads 
to the result already obtained in Ref. by a very different derivation. It should be remarked that our calculation 
shows that the HF approximation, which is equivalent to the so-called 'mean-field' description (see, e.g., 0), already 
includes nonperturbative effects. 

In general, the HF approximation for a contact interaction consists of writing the proper self energy as 



£(k,w„;^,T) = -^7^X3^ G(k',uv;/i,T) 



(31) 



with g = Anh 2 a/m |26|. Inserting this expression into Dyson's equation for the Green function and iterating the 
procedure we see that the Hartree-Fock approximation is equivalent to summing over all pure-bubble diagrams as 
mentioned above. Alternatively, we shall work out the bubble series directly. The starting point is the mean-field 
equation for the homogeneous density 

A^n hom = F 3/2 (-/3/i + 2/3gn hom ) (32) 
where F u (x) — (u(e~ x ), which is equivalent to the HF theory. With our scaled quantities, this equation reads 

f HF (A + v) = F 3/2 (A + v-f3 f i c + Aaf HF ) = F 3/2 (a + v + Aa[f HF (A + v) - C(3/2)]) (33) 



where we have inserted the lowest-order equation f3fi c = 4a((3/2) valid for small a |6j. The right-hand side is now 
expanded in powers of a: 



00 Li)" 

f HF (A + v)=J2 { —l-F 3/2 _ n (A + v)(4a) n [f HF (A + v)- C(3/2)] n . 



(34) 



Solving this implicit equation by iteration yields 

f H F(A + v) = F 3/2 (A + v)-4aF 1/2 (A + v)F 3/2 (A + v) 

F_ 1/2 (A + v)F 3 2 /2 (A + 



-(4a) 2 
-(4a) 3 



F 1 2 /2 (A + v)F 3/2 (A + v) 

-F? /2 (A + v)F 3/2 (A + v)- 
F_ 1/2 (A + v)F3 /2 (A + v) 



3 F_ 1/2 (A + v) F 1/2 {A + v) F 3 2 /2 (A + v) 



■0{a A 



(35) 
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where the subtracted function F v (x) = F v {x) — C(^) vanishes at x — for v > 1. The individual terms in i|35|) corre- 
spond to the HF-diagrams in the perturbation expansion for the homogeneous density (with the nonzero Matsubara 
frequencies taken into account). The first two terms once more give the previous expansion 1)30(1 . Using the Robinson 
expansion 

C„(e-*) = r(l - y)x v ~ x + ]T k), (36) 

k=0 



we find that the terms of order n > 2 in (|35[l diverge like a(a/yA + in the limit (A + v) — > 0, due to the leading 

Robinson term T(l — v)x" . Focusing attention upon this leading divergence, we replace each function F v (x) by its 
Robinson term. This step corresponds to dropping the nonzero Matsubara frequency contributions from the diagrams 
of Fig. |3 With the help of a computer algebra program, we obtain in this way from Eq. I|35|l the following series 
expansion of the divergent part of the density 

,(div), A , , _t^W , J[^W_ ^ 7/2 (4«) 6 57^(4^ 

Jhf 1 ' VAT^ 4(A + !I ) 3 / 2 8(A + vf/ 2 64(A + v)V 2 + ' ' ' ' 1 j 

Inserting this expression into Eq. (|24|) . we can easily perform the integration over v with the help of dimensional 
regularization. Note that it is crucial to carry out this ste p f or a finite A > 0, since otherwise all integrals would 
vanish according to Veltman's rule J dee a = for all a [Icl[l3j. The result of this calculation is 

»A) = r ( „ - i/2)r(-„)A" £ (r/ ' n + 1 ir (?? ' 1)?? (^f)" , (38) 

n>l ' ^ ' 

which is easily recognized as the series expansion of Y[r\ — l/2)r(— ry)[(167ra 2 + A) v — A 71 ]. In the critical limit A — > 0, 
this becomes 

I$Zk0) = (167r)T(„ - l/2)r(-r ? )a 2 ". (39) 

We see that due to the resummation procedure all diagrams in the perturbation series effectively contribute to the 
leading-order nonperturbative shift of the critical particle number. 

The same result l|39f) was also found in using a completely different approach. There, it was pointed out that 
this contribution can be derived from the behavior of the critical trapped density within a region around the trap 
center where < f3V(x) < a 2 . This is the nonperturbative regime by the Ginzburg criterion 20]. Here, in contrast, 
we use the perturbation expansion, which is valid far away from the trap center, to obtain the same result. 

We now turn to evaluating N ^ (T, /x) , i.e., the second contribution to the trapped particle number which is deter- 
mined by the diagrams displayed in Fig.0 For the calculation, we shall make use of the high-order perturbative loop 
expansions that were derived in Refs. [lall7| . These expansions allow us to evaluate the most divergent contributions 
n hom t° tne homogeneous density, as defined above. Proceeding along the lines indicated in the study of 
we use fi^om to first obtain a series expansion for 7V( 2 ' dlv ) and then to resum the series to find the change in the 
critical particle number. This time the resummation cannot be performed exactly, which leads us to apply variational 
perturbation theory (VPT). 

A problem in this procedure concerns the fact that Refs. do not calculate nj^'^ in terms of /i as we need it, 

but rather as a function of £ = — (3[fi — ?i£(0, fi)] (in other words, we would need the diagrammatic expansion leading 
to Fig. |3 to be carried out similar to Fig- El i- e -? including bubble contributions). In the following we will ignore this 
difference and approximate the exact expression 

/>OC 

N (2,dW (A = pfa _ m)) = CQnst . x / d£ £l) -3/2 A | n (2,div) ^ _ A _ p£) (4Q) 

Jo 

by 

tf(W (A )« const, x / de'^ 3/2 A^of V) (A + r). (41) 
Jo 

In this equation n^'^ denotes the divergent homogeneous density as a function of £, i.e., ^om^ (/3/z) = 



"horn 



(2 ' div) (£(/?M))- In Refs. E1IT3, n hom V ^ (£) * s calculated in terms of a high-order perturbative loop expansion: 

4*2^(0 = aX>(^) • (42) 
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FIG. 4: Results for C^irf) from VPT: (bold curve) fixed uo — 0.805; (dash-dotted) self-consistent determination of oj' . Note the 
divergence as r) — * 1. 



From power counting, we expect /5?il](0, A) to have an expansion of the form a 2 si (a/VA)' 2 for the leading-order 
divergence. Using this expansion we see that the exact expression for iV^ 2 ) and our approximation differ somewhat in 
the higher-order coefficients of the weak-coupling expansion. Since the results of the resummation are most strongly 
affected by the lower-order coefficients, we neglect the error introduced in this way. We expect that this simplification 
will not change the main features of our results. 

We thus insert the expansion Q42f) into Eq. (|41|l . The coefficients bi appearing in (|42|) are related to the a/'s of Ref. 
[ijj by bi = 3847r 3 (247r)'- 2 a ; . Applying dimensional regularization, we now perform the integration 



V 



i-2 



a ,fj±y"«H (43) 



This result constitutes a weak-coupling expansion of 7V( 2 ' dlv )(T, /i). However, we have to consider the limit A 
which is a strong-coupling limit. We assume that the integrand in l|43(l has a strong-coupling expansion [Tfl 



rfa-i/2) 



I dvv^A* T n^\A + v)=a^Y,B m ^ . (44) 



The leading power of a is fixed by dimensional considerations, i.e., by formally equating (|43|l with (J33J, since the 
result can only depend on the parameter a/\f~K [also see the discussion following Eq. Q14[l] . The subleading powers 
are multiples of the universal Wegner exponent governing the approach to scaling as shown in |l6l Il7| . This exponent 
reflects the influence from the anomalous dimensions of quantum field theory. 

Our task is now to determine the coefficient Bq of the strong-coupling expansion. This will yield the leading-order 
contribution iV^ 2,dlv ^ to the shift in N. The result is obtained from field-theoretic VPT in analogy to the calculations 
for the homogeneous system in Ref. |l6j| . As in that reference we perform two variants of the calculation, one by fixing 
lj' to have the known value 0.805, and one by determining w' order by order self-consistently. At first glance, this 
approach still leaves grounds for scepticism, since the perturbation expansion obtained from the LDA is divergent 
in two ways. First, the expansion coefficients grow factorially, so that the radius of convergence is zero. Second, 
the series has to be evaluated at a very large argument which goes to infinity when approaching the critical point. 
Fortunately, these two unpleasant properties are familiar from perturbation expansions of critical exponents, for which 
it has been shown that a resummation by field-theoretic VPT leads to the correct results. We can thus also safely use 
this method here. 

Within our approach, the change in the critical particle number is finally obtained as an expansion 

1 I v (Q) = ari + l) + C 1 {r l )*+[C'M+C' 2 {r l )\a 2v + G(a 2 ) (45) 



rfa-1/2) 



with Ci(r]) determined from Eq. (|30|) and C[(rj) — (167r)T(— rj) as in Eq. (|39|) . The coefficient C' 2 (ri) emerges from 
VPT. In Fig. 21 we show the VPT calculation for C' 2 (rf) using the seven-loop data from Ref. [17|. The result with 
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FIG. 5: Shift of critical temperature as a function of the potential shape parameter r] for fixed a/Ar = 10~ 4 : (full curve) Result 
including all terms in 14511 and 62(17) calculated with fixed u' = 0.805; (dash-dotted) same as full curve but with u>' determined 
self-consistently; (dashed) 62(77) omitted (i.e., "mean-field result"); (dotted) C[(rj) and 62(77) omitted (i.e., only perturbative 
term) . 



self-consistent determination of ui' (bold curve) is compared to the calculation with a fixed value of ui' — 0.805. We see 
that both curves agree reasonably well; the calculation of C" 2 (r)) thus does not depend too sensitively on the exponent. 
For ?7>0.94, the self-consistent calculation of lo' does not converge anymore; one would probably have to extend the 
resummation to higher orders to resolve this issue. More importantly, however, we see that the results for C' 2 (rj) 
diverge in the limit of 77 — > 1. This behavior can be traced back to the presence of the T functions in Eq. Ij43(l . The 
same divergence also appears in the coefficient C[ (77) . We will discuss the significance of these divergences below in 
connection with the calculation of the critical temperature. 

In Fig. we plot the final result, i.e., the shift in the critical temperature, as determined from Eqs. (|25|l and Q45p. 
The figure reflects the characteristic qualitative features discussed in Sec. Ill and confirms our previous conclusions. 
The shift is displayed as a function of 77 for a fixed value of a = a/\r — 10 -4 . The full and the dash-dotted curves 
show the complete result including all terms from Eq. (|45H . The two curves are respectively based on the calculation 
of C2 (77) with u' kept fixed or determined self-consistently. The discrepancy between these results can be considered 
an estimate for the error in the calculation of T c that is introduced by the resummation procedure. Since the curves 
are almost indistinguishable, we can conclude that the resummation contributes only a small inaccuracy in addition 
to potential further errors introduced by the other approximations. 

For the dashed curve in Fig. the term C 2 (r])a 2v has been omitted, it thus displays the mean-field result below 
second order 0. The dotted curve shows the perturbative contribution due to Ci{rf)a. We see that the full and 
the mean-field result closely approach the perturbative first-order approximation long before values of 77 around 1 
are reached. The divergence around 77 = 1 introduced by the behavior of C[ and C 2 is restricted to a very small 
interval. Thus it is reasonable to expect that the true behavior of T c in this small regime remains well described by 
the first-order approximation and that the divergence is only an artifact of our calculation. Furthermore, in Ref. it 
was shown that within mean-field theory this divergence is compensated for by a pole in the second-order contribution 
to the critical particle number. We can also expect the same behavior in the present case, i.e., beyond mean-field 
theory. 



V. SUMMARY AND CONCLUSIONS 



We have calculated the shift of the critical temperature of interacting Bose gases trapped in a general power- 
law potential of the type Ui\xi\ Pi . The objective was to understand how this shift changes when we pass from 
homogeneous to harmonically trapped systems, interpolating between these limits by changing the power of the 
potential. While the homogeneous case is influenced strongly by nonperturbative critical fluctuations, the harmonic 
case can be calculated pcrturbatively. 

We have restricted our attention to the thermodynamic limit, which allowed us to use the local-density approxima- 
tion in which the Bose gas is assumed to be locally homogeneous at each point. In Sec. Ill, we gave justification for 
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this procedure. The main result is the shift formula 




+ 0(a 2 ). 



(46) 



It contains a linear, perturbative part, which is relevant for all potentials, and a nonperturbative contribution pro- 
portional to (a/ At) 211 ■ For small a/ At, the latter contributes significantly only in the quasihomogeneous regime. 

The presence of the (a/ At) 211 term was derived from scaling and resummation arguments, and we gave a simple 
physical explanation for its appearance based on the Ginzburg criterion pol ]. Our results show how the growing 
inhomogeneity of the potential reduces the significance of critical fluctuations. 

We also performed an explicit calculation of the nonperturbative coefficient D'(rj). In spite of the approximate 
character of the calculation, the result should be quite accurate. Our approach was based on the resummation of 
divergent perturbation series with the help of field-theoretic variational perturbation theory. In the course of the 
derivation, it was shown that the simple Hartree-Fock approximation also contains a nonperturbative contribution. 

Finally, our study indicates that the higher-order T c shift in traps with r) = 1, for instance in J^. l/|a;i| 6 potentials, 
should be particularly interesting and difficult to investigate. In the limit r\ — > 1, we find a divergence in the coefficient 
D'(rj) which governs the nonperturbative contribution proportional to (a/ At) 2 - It remains to be investigated whether 
this divergence is canceled by other genuinely second-order terms, as it is expected from mean-field theory. For 
future work, it might also be of interest to study the influence of other types of external potentials on the transition 
temperature, such as optical lattices |28j. 
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